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Abstract 

qq . We study the behavior of simple principal pivoting methods for the P-matrix 

linear complementarity problem (P-LCP). We solve an open problem of Morris by 
showing that Murty's least-index pivot rule (under any fixed index order) leads to a 
quadratic number of iterations on Morris's highly cyclic P-LCP examples. We then 
show that on K- matrix LCP instances, all pivot rules require only a linear number of 
iterations. As the main tool, we employ unique-sink orientations of cubes, a useful 
C"| , combinatorial abstraction of the P-LCP. 

1 Introduction 

The third author of this paper still vividly recollects his visit to Victor Klee at the 
*^ . University of Washington (Seattle) in August 2000. 

£\i \ We had a memorable drive in Vic's car from the hotel to the math depart- 

ment. Vic skipped all the small talk and immediately asked: "Do you think 

t — that the simplex method is polynomial?" I still remember how awkward I 

felt to be asked this question by the very person who had provided the first 
and seminal insights into it. My (then as now quite irrelevant) opinion shall 
remain between Vic and me forever. 

^ The seminal work referred to above is of course the 1972 paper of Klee and Minty 

entitled "How good is the simplex algorithm?" [TU] . It deals with the number of itera- 
tions that may be required in the worst case to solve a linear program (LP) by Dantzig's 
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simplex method [9], which was at that time the only available practical method. Klee 
and Minty exhibited a family of LPs (now known as the Klee-Minty cubes) for which 
the number of iterations is exponential in the number of variables and constraints. 

Dantzig's specific pivot rule for selecting the next step is just one of many conceivable 
rules. The question left open by the work of Klee and Minty is whether there is some 
other rule that provably leads to a small number of iterations. This question is as open 
today as it was in 1972. 

Nowadays, the simplex method is no longer the only available method to solve LPs. In 
particular, there are proven polynomial-time algorithms for solving LPs — Khachiyan's 
ellipsoid method [18J and Karmarkar's interior point method [T7] — that are based on 
techniques developed originally for nonlinear optimization. It is still unknown, though, 
whether there is a strongly polynomial-time algorithm for solving LPs, that is, an algo- 
rithm for which the number of arithmetic operations does not depend on the bit lengths 
of the input numbers. 

The P-matrix linear complementarity problem. In this paper we are concerned with 
pivoting methods for the P-matrix linear complementarity problem (P-LCP), a promi- 
nent problem for which neither polynomial-time algorithms nor hardness results are 
available. 

In general, an LCP is given by a matrix M € M nxn and a vector q € K n , and the 
problem is to find vectors w, z € R n such that 

w-Mz = q, w,z>0, w T z = 0. (1) 

It is NP-complete to decide whether such vectors exist [TJ [20] . However, in a P-LCP, 
M is a P-matrix (meaning that all principal minors are positive), and then there are 
unique solution vectors w,z for every right-hand side q [30]. The problem of efficiently 
finding them is unsolved, though. 

Megiddo has shown that the P-LCP is unlikely to be NP-hard. For this, he considers 
the following (more general) variant: given M and q, either find w,z that satisfy (PQ), 
or exhibit a nonpositive principal minor of M. NP-hardness of the latter variant would 
imply that NP = coNP [23]. 

There is a different notion of hardness that might apply to the P-LCP as a member of 
the complexity class PPAD [28]. This class has complete problems, and no polynomial- 
time algorithm is known for any of them. It has recently been shown that the famous 
problem of finding a Nash equilibrium in a bimatrix game is PPAD-complete [6]. Some 
researchers consider this to be an explanation why (despite many efforts) no polynomial- 
time algorithm has been found so far. Incidentally, this is the second problem for which 
Megiddo showed in his technical report [23] that it is unlikely to be NP-hard. However, 
it is currently not known whether the P-LCP is also PPAD-complete. 

There are various algorithms for solving P-LCPs, among them the classic method 
by Lemke [22] as well as interior point approaches |214 I2U| . Still, there is no known 
polynomial-time algorithm for P-LCPs. For example, the complexity of interior point 
algorithms depends linearly on a matrix parameter k [20] that is not bounded for P- 
matrices. 



In this work, we focus on combinatorial methods for P-LCPs and in particular on 
simple principal pivoting methods that share their essential idea with the simplex method. 

Simple principal pivoting methods. Let B C {1,2, .. . , n}, and let Ab be the n x n 

matrix whose ith column is the ith column of — M if i £ B, and the zth column of the 
n x n identity matrix I n otherwise. If M is a P-matrix, then Ab is invertible for every 
set B. We call B a basis. If A^ q > 0, we have discovered the solution: let 

jo ifieB j(A^q)i iiieB 

w« := < , , Zi := < • (2) 

\{A~ l q)i \iiiB [0 iliiB 

If on the other hand Ag q ^ 0, then w and z defined above satisfy w — Mz = q and 
w z = 0, but w, z > fails. In principal pivoting, one tries to improve the situation by 
replacing the basis B with the symmetric difference B © C, where C is some nonempty 
subset of the "bad indices" {i : (Ag q)i < 0}. The greedy choice is to let C be the set 
of all bad indices in every step. For some P-LCPs, however, this algorithm cycles and 
never terminates [27 |H In simple principal pivoting, C is of size one, and a pivot rule 
is employed to select the bad index. Simple principal pivoting methods are sometimes 
called Bard-type methods and were first studied by Zoutendijk [31] and Bard [2]. For a 
detailed exposition on simple principal pivoting methods see, for instance, [HI Section 4.2] 
or [271 Chapter 4]. 

The digraph model and unique sink orientations. There is a natural digraph model 
behind principal pivoting, first studied by Stickney and Watson [31J . The graph's vertices 
are all bases B, with a directed edge going from B to B' HB@B' = {i} and (A] 3 1 q)i < 0o 
The underlying undirected graph is the graph of the n-dimensional cube. A principal 
pivot step considers the cube that is spanned by the outgoing edges at the current vertex, 
and it jumps to the antipodal vertex of some subcube of it. A simple principal pivot 
step can be interpreted as the traversal of an outgoing edge, and this is the analogy with 
the simplex method for linear programming. 

The above digraph has the following specific property: every nonempty subcube (in- 
cluding the whole cube) has a unique sink, and the global sink corresponds to the solution 
of the P-LCP [3T]. In the terminology of Szabo and Welzl [32], we are dealing with a 
unique-sink orientation (USO). 

The goal of this paper is to deepen the understanding of the USO description of simple 
principal pivoting methods for the P-LCP. On a general level, we want to understand 
whether and how algebraic properties of a P-LCP (which are well studied) translate to 
combinatorial properties of the corresponding USO (which are much less studied), and 
what the algorithmic implications of such translations are. The principal motivation 



x For n — 3, one such example is the P-LCP (JU whose digraph model (see next paragraph) appears in 
Figure [T] Cycling occurs for any of the three start bases just below the top vertex in Figure [T] 

2 We require q to be nondegenerate, meaning that (Ag 1 q)i ^ for all B and i. This is no serious 
restriction since a nondegenerate problem can always be obtained from a degenerate one by a symbolic 
perturbation of q. 



behind this research is the continuing quest for a polynomial-time P-LCP algorithm. In 
this paper, the concept of a USO serves as the main tool to obtain two more positive 
results for simple principal pivoting. 

We believe that this combinatorial approach has some untapped potential for the 
theory of (simple) principal pivoting methods for LCPs. 

The (randomized) method of Murty and the Morris orientations. The simple prin- 
cipal pivoting method of Murty (also known as the least-index rule) [25] works for all 
P-LCP instances, but it may take exponentially many iterations in the worst case [26] . 
As a possible remedy, researchers have considered randomized methods. The two most 
prominent ones are RandomizedMurty (which is just the least-index rule, after ran- 
domly reshuffling the indices at the beginning), and RandomEdge (which performs a 
purely random walk in the USO). However, Morris found a family of P-LCP instances (we 
call their underlying digraphs the Morris orientations) on which RandomEdge spends 
a long time running in cycles and thus performs much worse than even the exhaustive 
search algorithm [24] . 

For RandomizedMurty there are as yet no such negative results on P-LCP in- 
stances. In particular, on Murty's worst-case example, this algorithm takes an expected 
linear number of steps [TT]; the expected number of steps becomes quadratic if we allow 
arbitrary start vertices [TJ. 

On general P-LCP instances, the expected number of iterations taken by Random- 
izedMurty is subexponential if the underlying digraph is acyclic [H], but it is unknown 
whether this also holds under the presence of directed cycles (as they appear in the 
Morris orientations, for example). 

We prove that not only the randomized variant, but actually any variant of the least- 
index rule obtained by some initial reshuffling of indices leads to a quadratic number 
of iterations on the Morris orientations. In particular, this "kills" another family of 
potentially bad instances for the RandomizedMurty rule. 

K-matrix linear complementarity problems. A K-matrix is a P-matrix so that all 
off-diagonal elements are non-positive. LCPs with K-matrices (K-LCPs) appear for 
example in free boundary problems [Q Section 5.1] and in American put option pricing 
[3j E]. It is known that every K-LCP instance can be solved in polynomial time (even 
by a simple principal pivoting method [E1E9], see also [HJ Section 4.7]), but we prove 
something stronger: every simple principal pivoting method takes only a linear number 
of iterations. We obtain this result by showing that in K-LCP-induced USOs, all directed 
paths are short. Our approach is to extract combinatorial structure from the algebraic 
properties of K-matrices. Subsequently, we use the distilled combinatorial information 
to get our structural results. 

LCPs with sufficient matrices and the criss-cross method. Let us step back and take a 
broader view that brings together LP, the P-LCP, pivoting methods, and computational 
complexity: linear complementarity problems with sufficient matrices generalize both 



LP and the P-LCP while still being amenable to pivoting approaches. The criss-cross 
method [12^113] may be considered an extension of Murty's algorithm. In fact, the results 
of |12j make LCPs with sufficient matrices the largest known class of LCPs for which 
Megiddo's techniques of [23] still apply, so that NP-hardness is unlikely. Determining 
the complexity of sufficient matrix LCP remains a tough challenge for the future. 

2 Unique-sink orientations 

Now we formally introduce the digraph model behind P-LCPs and show some of its basic 
properties. The model was first described by Stickney and Watson in 1978 [31] . 

We use the following notation. Let [n] := {1, 2, . . . , n}. For a bit vector v & {0, l} n 
and I C [n], let v © / be the element of {0, l} n defined by 

This means that v © I is obtained from v by flipping all entries with coordinates in /. 
Instead of v © {i} we simply write v © i. 

Under this notation, the (undirected) n-cube is the graph G = (V, E) with 

V :={0,l} n , E:={{v,vQi} :v£V, i€ [n]}. 

A subcube of G is a subgraph G' = (V , E 1 ) of G where V = {v © / : / C C} for some 
vertex v and some set C Q [n] , and E' = E n ( 2 ) . 

2.1 Definition. Let (j) be an orientation of the n-cube (a digraph with underlying 
undirected graph G). We call <p a unique-sink orientation (USO) if every nonempty 
subcube has a unique sink in (p. 

Figure Q] depicts a USO of the 3-cube. The conditions in Definition 12.11 require the 
orientation to have a unique global sink, but in addition, all proper nonempty subcubes 
must have unique sinks as well. In Figure Q] there are six 2-dimensional subcubes, twelve 
1-dimensional subcubes (edges) and eight 0-dimensional subcubes (vertices). For edges 
and vertices, the unique-sink condition is trivial. The figure also shows that USOs may 
contain directed cycles. 

If V is the vertex set of a subcube (as defined above), then the directed subgraph of <f) 
induced by V is denoted by 0[V']. 

It is convenient to view the orientation of the n-cube G as a mapping <j) : {0, l} n — > 
{— , +} n , where (fi(v)i = — if the edge between v and v © i is oriented towards v © i, and 
4>(v)i = + if it is oriented towards v. To encode a vertex v along with the orientations of 
its incident edges, we can then use a configuration of n bits and n signs, where the ith 
sign is 4>{v)i (see Figured]). If an orientation 4> of G contains the directed edge [v , v ©i), 

we write v — > v © i, or simply v — > v © i if 4> is clear from the context. 
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Figure 1: A USO of the 3-cube; for each vertex, the orientation of incident edges is 
encoded with a sign vector in { — ,+}" (— for an outgoing edge; + for an 
incoming edge). The unique global sink and a directed cycle are highlighted. 

Let (f) be an orientation of the ra-cube and let F C [n]. Then cf>( ' is the orientation of 
the n-cube obtained by reversing all edges in coordinates contained in F; formally 

6(F) I v -^> v i if i & F, 

v > v © i :<3> < , 

\v®i-*v if i G F. 

2.2 Proposition ([32]). If (j) is a USO, then <f)( F > is also a USO for an arbitrary sub- 
set FC[n]. D 

This proposition can be proved by induction on the size of F; it essentially suffices to 
show that under cp^ 1 ' > , the whole cube still has a unique sink. 

Now we can derive the unique- completion property. Its meaning is that we can pre- 
scribe the bits for some coordinates and the signs for the remaining coordinates, and in 
a USO there will always be a unique vertex that satisfies the prescription. In particular, 
it follows that any USO (and any of its subcubes) also has a unique source. 

In fact, the unique-completion property characterizes unique-sink orientations. 

2.3 Lemma. An orientation 4> of the n-cube G is a unique-sink orientation if and only 
if for every partition [n] = AuB and every pair of maps a : A — > {0, 1} ; : B — » {— , +} 



there exists a unique vertex v such that 

Vi = a(i) for all i G A, (3a) 

( j ) ( v ) i = p(i) for all i G B. (3b) 

Proof. First suppose that is a USO. Let U := {u : Ui = a(i) for all i G ^4} and consider 
the subcube <j)\U\. Let F := {i G 1? : /3(i) = — }. As we have just noted, 4>( F ' is a unique- 
sink orientation, and thus the subcube 4>[U] has a unique sink v with respect to <jy- F h 
This vertex u is the only vertex that satisfies (|3H) and (J3b). 

Conversely, let be an orientation of the n-cube which satisfies the unique-completion 
property and let U := {u ® I : I C C}. Set B := C and A := [n] \ C, and let a(i) := ui 
for i G A and /3(i) := + for i & B. Then the unique vertex v satisfying ([3^,) and §3]p) is 
the unique sink of the subcube <f>[U]. D 

P-LCP-induced USOs. We now formally define the USO that is induced by a nondegen- 
erate P-LCP instance P-LCP(M,<?) [31]. For v G {0, l} n , let B(v) := {j G [n] : v 5 = 1} 
be the canonical "set representation" of v. Then the unique-sink orientation eft induced 
by P-LCP (M,q) is 

v^v®i :& (A-^ v) q)i<0. (4) 

Recall that Aq is the (invertible) n x n matrix whose ith column is the ith column 
of — M if i G B, and the ith column of the n x n identity matrix I„ otherwise. 

In this way, we have reduced the P-LCP to the problem of finding the sink in an 
implicitly given USO. Access to the USO is gained through a vertex evaluation oracle 
that for a given vertex v G {0, 1}™ returns the orientations (f)(v) G {— , +} n of all incident 
edges. In this setting, the aim is to find the sink of the USO so that the number of 
queries to the oracle would be bounded by a polynomial in the dimension n. In the case 
of a P-LCP-induced USO, a (strongly) polynomial-time implementation of the vertex 
evaluation oracle is immediate from §%§ . From the sink of (j), we can reconstruct solution 
vectors w and z as in ^). The fact that the n-cube orientation defined in this way is 
indeed a unique-sink orientation was proved by Stickney and Watson [31J. 

Not every USO is P-LCP-induced; for n = 3, the P-LCP-induced USOs are charac- 
terized in [31] . but for n > 4, we have only a necessary condition. This condition was 
originally proved by Holt and Klee for polytope graphs oriented by linear functions |16j . 

2.4 Theorem f[15|). Every P-LCP-induced USO of the n-cube satisfies the Holt-Klee 
condition, meaning that there are n directed paths from the source to the sink with pair- 
wise disjoint sets of interior vertices. □ 

3 Pivot Rules 

In the USO setting resulting from P-LCPs, simple principal pivoting can be interpreted 
as follows: start from any vertex, and then proceed along outgoing edges until the global 
sink (and thus the solution to the P-LCP) is reached. A pivot rule R determines which 



outgoing edge to choose if there is more than one option. Here is the generic algorithm, 
parameterized with the rule R. It outputs the unique sink of the given USO <f>. 

3.1 Algorithm. 

SimplePrincipalPivoting r (0) 
Choose v G {0, l} n 

while O := {j G [n] : v -> v © j} / do 

choose i G O according to the pivot rule R 

v := v © i 
end while 
return v 

Note that when eft contains directed cycles as in Figure [TJ this algorithm may enter 
an infinite loop for certain rules R, but we consider only rules for which this does not 
happen. In the USO setting we are restricted to combinatorial rules. These are rules that 
may access only the orientations of edges {v, v © j} (as given by the vertex evaluation 
oracle) but not any numerical values such as {A~Z, ^q)j that define these orientations. 

Let us now introduce the combinatorial pivot rules that are of interest to us. We 
write them as functions of the set O of candidate indices. The first one is Murty's least- 
index rule Murty that simply chooses the smallest candidate. In the cube, this may be 
interpreted as a "greedy" approach that always traverses the first outgoing edge that it 
finds (in the order of the cube dimensions). 

Murty(O) MURTY 7r (C) 

return min(O) return 7r(iniri(7r -1 ((9))) 

It is easy to prove by induction on n that this rule leads to a finite number of iterations, 
even if the USO contains directed cycles. The rule has a straightforward generalization, 
using a fixed permutation n G S n : choose the index ir(i) with the smallest i such that 
v — > v © 7r(i). This just reshuffles the cube dimensions by applying n, and with ir = id, 
we recover Murty's least-index rule. 

The natural randomized variant of Murty,,- is RandomizedMurty: at the begin- 
ning, choose the permutation it uniformly at random from the set S n of all permutations, 
and then use Murty^- throughout. 

Finally, RandomEdge stands for the plain random walk: in each iteration, simply 
choose a random candidate. 

RandomizedMurty(O) RandomEdge(O) 

if called for the first time then choose i G O uniformly at random 

choose 7r G S n uniformly at random return i 

end if 
return Murty^-(O) 

Unlike the previous rules, RandomEdge may lead to cycling in SimplePrincipal- 
PivotinGr, but the algorithm still terminates with probability 1, and with an expected 



number of at most (n + l)n n iterations. This is a simple consequence of the fact that 
there is always a short directed path to the sink. 

3.2 Lemma ([31, Proposition 5]). Let t E {0, l} n be the global sink of an n-cube USO 4>, 
and let v E {0, l} n be any vertex. If k is the Hamming distance between v and t, then 
cj) contains a directed path of length k from v to t. □ 

4 The Morris orientations 

Morris proved that under R = RandomEdge, Algorithm 13.11 may be forced to perform 
an expected exponential number of iterations [24J. More precisely, at least ((n — l)/2)! 
iterations are required on average to find the sink of the n-cube USO generated (as 
described in Sectional) by the LCP(M, q), with 
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Here, n must be an odd integer in order for M to be a P-matrix. 

To prove this result Morris first extracted the relevant structure of the underlying 
USO, and then showed that on this USO, RandomEdge runs in cycles for a long time 
before it finally reaches the global sink. It was left as an open problem to determine the 
expected performance of RandomizedMurty, the other natural randomized rule, on 
this particular USO. 

We solve this open problem by showing that for any permutation n, Murty^ incurs 
only 0(n 2 ) iterations on the n-dimensional Morris orientation. In particular, this bound 
then also holds for RandomizedMurty. 

4.1 Theorem. For any ir € S n , Alaorithm \3.1\ with R = MuRTYjr finds the sink of the 
n-dimensional Morris orientation after at most 2ra 2 — (5ra — 3)/2 iterations, from any 
start vertex. 

Computational experiments suggest that the worst running time is attained for the 
identity permutation if the algorithm starts from (the global source). Therefore we 
analyze this case thoroughly and compute the precise number of iterations. 

4.2 Theorem. Algorithm \3.1\ with R = MuRTYij finds the sink of the n-dimensional 
Morris orientation after (n 2 + l)/2 iterations, starting from 0. 

Experimental data suggest that this is an upper bound on the number of iterations 
for any permutation -k and start vertex 0. Allowing start vertices different from 0, we 



get slightly higher iteration numbers, but n = id continues to be the worst permutation. 
In view of this, we conjecture that n 2 /2 + 0(n) iterations actually suffice for all tt and 
all start vertices. This would mean that the bound in Theorem 14. II is still off by a factor 
of 4. 

The combinatorial structure. The n-cube USO resulting from ([5]) can be described 
in purely combinatorial terms |24} Lemma 1]. Here we make its structure even more 
transparent by exhibiting a finite state transducer (finite automaton with output) with 
just two states that can be used to describe the USO, see Figure EJ For the remainder 
of this section, we fix n to be an odd integer and we are considering the n-dimensional 
Morris orientation 6. 




Figure 2: The finite state transducer that generates the Morris orientations. Each tran- 
sition is labeled with an input symbol (bit) and an output symbol (sign) 

The orientation is then determined for each vertex v = (v\,V2, ■ ■ ■ ,v n ) G {0, l} n as 
follows. If v\ = i)2 = ■ ■ ■ = v n = 1, then v is the global sink. Otherwise choose 
some i in [n] so that V{ = and consider the bit string t>j_i?jj_2 . . . viv n v n -\ . . . vi as 
the input string to the transducer. In other words, the transducer is reading the in- 
put from right to left, starting immediately to the left of some 0. The output string 
4>{v)i-i4>{v)i-2 ■ ■ ■ 4>(v)i4>(v) n (p(v) n -i . . . <f>(v)i then determines the orientation of each 
edge incident to v. The choice of i does not matter, because after reading any zero, the 
transducer is in the state S. 

For instance, let n = 5 and let v = (1,0, 1, 1,0). Then the transducer takes v^ = 1, 

outputs + and switches to state T; reads v% = 1, outputs — and switches to S; reads 

V2 = 0, outputs — and stays in S; reads v\ = 1, outputs + and switches to T; and 

finally reads v§ = 0, outputs + and switches to S. The resulting configuration for v is 
( 1 o 1 l o \ 

Figure Q] depicts the Morris orientation for n = 3. It can easily be checked that 
the above procedural definition of the orientation is equivalent to the one given in |24[ 
Lemma 1]. In particular, we obtain a USO. It has the remarkable symmetry that the 
set of configurations is closed under cyclic shifts. 

Pivoting. From the transducer we can easily derive the sign changes in the configura- 
tion ( 0^) ) that are caused by a pivot step v — > u := v © {i}. We always have <j){v)i = — 
and 4>{u)i = +. 
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If Vi = 1, we get Ui = 0. In processing v and u, the transducer then performs exactly 
the same steps, except that at the ith coordinate different transitions are used to get 
from T to S. But this implies (j>(u)j = <f>(v)j for j ^ i. Here is an example for such a 
pivot step (with i = 3; affected signs are indicated): 

(i£U*)-(i£gU) (6) 

If Vi = 0, we get Ui = 1, meaning that at the ith coordinate the transducer stays in S 
for v but switches to T for u (assuming that u is not the sink). This implies that signs 
change at all coordinates i — 1, i — 2 down to (and including) the next coordinate k (wrap 
around possible) for which V}~ = u& = 0. In our example we have such a step for i = 2: 

[HI 1 B\ ,/l 11 10\ (7\ 

\+ h + ) ~~ * \e e - + e) \' ) 

In both v and u, the signs in the block of l's at indices k + 1, . . . , i — 1 alternate. 

Levels. It will be useful to define the following function on the vertices of the cube. 
Let £{v) be the number of coordinates where (_) appears. Formally, 

£(v) := \{i : Vi = and (/>(«)* = -} | . 

The number £(v) is called the level of v. From the two types of pivots, it is easy to see 
that the value of £ does not increase along any directed edge [33, Lemma 3]. Indeed, if 
u is an out-neighbor of v, then either £(u) = £{v) (this happens in every pivot of type 
([6]), and in pivots of type (JTj) with an odd block of l's at indices k + 1,. . . ,i — 1), or 
£ decreases to the next possible value: £(u) = £{v) — 2, or £{u) = if £(v) = 1. In 
particular, the values of £ lie in the set {n, n — 2, ... ,3,1, 0}. The sink is the only vertex 
at level 0. 

Let us make a small digression and briefly explain why RandomEdge is slow on (j>. 
Let L(v) := \{i : Vi = 1 and 4>(v)i = — }|. So the outdegree of v is £{v) + L(v). Now let 
v be a vertex at level 1. In pivots of type ©, L decreases by one, whereas in pivots of 
type ([7]) that do not reach the sink, L increases by one. The latter occurs if and only if 
RandomEdge pivots at the unique (_), which happens with probability l/(L(v) + 1). 

At level 1, we can thus interpret RandomEdge as a random walk on the integers, 
where k — » k + 1 with probability l/(k + 1), and k — > k — 1 with probability k/(k + 1). 
Therefore the walk is strongly biased towards 0, but in order to reach a neighbor u of 
the sink, we need to get to k = L(u) = (n — l)/2. This takes exponentially long in 
expectation 



A quadratic bound for Murty,,.. To prove Theorem 14.11 we first derive a (somewhat 
surprising) more general result: for any pivot rule, starting from any vertex v with 
Vk = for some k, it takes 0(n 2 ) iterations to reach a vertex u with uj- = 1. This can be 
viewed as a statement about the USO induced by the n-dimensional Morris orientation 
on the cube facet {v : v^ = 0}. The statement is that in this induced USO, all directed 
paths are short (in particular, the induced USO is acyclic). 
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4.3 Lemma. Let v € {0, l} n be a vertex, and let k be any index for which Vk = 0. Start- 
ing from v, Algorithm ] 3. 1\ (with an arbitrary pivot rule r) reaches a vertex u satisfying 
Ufc = 1 after at most ( g ) iterations. 

Proof. We employ a nonnegative potential that decreases with every pivot step. By a 
cyclic shift of coordinates we may assume that k = 1, and thus v\ = 0. 
Let us define 

Nq(v) := {j € [n] ■v j = and <f>(v)j = -}, 
JVi(«) := {j G [n] : uj = 1 and 0(f) j = -}. 

The potential of t; is the number 

p(«) = |JNTi(v)|+ Y, J- 

j£N {v) 

Now we will observe how the potential changes during a pivot step v — > u = v © {i}. 
We always have <^(i>)j = — , </>(«)« = +• There are two cases. If v% = 1, as in ([6]), the 
configurations of v and u differ only at index i. Hence Nq(u) = Nq(v) and N\(u) = 
Ni(v) \ {i}. Therefore p(u) = p{v) — 1. 

In the other case we have Vi = 0, as in ([7]). If % = 1, we are done, since then ui = 1. 
Otherwise, v is of the form 

V = (VI,..., Vj , Vj + l,. . .,Vj-l , Vi ,v i+ i,...,v n ), 

* w^ ° 

with 1 < j < i. If i — j is odd, there is an even number of l's between Vj and Vi, implying 
that |iVi(«)| = \Ni(v)\ and N (u) = N (v) \ {j,i}. We then have p{u) < p(v). lii—j is 
even, there is an odd number of l's between Vj and Vi, implying that |iVi(u)| = \Ni(v)\ + l 
and Nq(u) = Nq(v) U {j} \ {i}. Since i — j > 2 in this case, we get p(u) < p(v) — 1. 

To summarize, the potential decreases by at least 1 in every pivot step v := v(Bi, and 
as long as v is not the sink, we have p(v) > 1. The highest potential is attained by the 
source, p(0) = ( n ^ ). Therefore after at most ( n ^" ) steps we have v\ = 1. □ 

If w happens to be at level 1 already, we obtain a better bound, since |iVo(v)| = 1 
in this case. In fact, the vertex with largest potential in level 1 such that v\ = is 
v = ( j i 1 i •■■ j i 1 ° ) , and its potential is p(v) = 3(n - l)/2. 

4.4 Corollary. Let v € {0, l} n be a vertex at level 1, and let k be any index for which 
Vk = 0. Starting from v, Algorithm \3.1\ (with arbitrary pivot rule r) reaches a vertex u 
satisfying u^ = 1 after at most 3(n — l)/2 iterations. □ 

Now MuRTY^ comes in: we will apply either Lemma 14.31 or Corollary 14.41 to vertices 
v with v v u\ = and ?V(j+i) = • • • = Uwra) = 1, for some i. The next lemma shows that 
the l's at coordinates Tv(i + 1), . . . , n(n) will stay. 

4.5 Lemma. Let i € [n]. After Algorithm ic!. 1\ with R = Murty^ visits a vertex v with 
v Tr(i+i) = • • • = %( n ) = 1; ^ wi// visit only vertices u satisfying u % u + u = • • • = u^( n \ = 1. 
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Proof. At every iteration, the algorithm replaces the current vertex v with v © tt(j), 
where j is the smallest index such that v —> v (B ir(j). It follows that coordinates 
n(i + 1), . . . , 7r(n) will be touched only if v © tt(j) —> v for all j < i. When this holds for 
the first time, we have 

v *(j) = 1 ' for 3 > *> 
0O)tt(.j) = +' for j < i. 

The unique-completion property (Lemma I2.3P then implies that we have already 
reached the sink. □ 

Now we can put it all together. 



Proof of Theorem \4-l\ Let v be the start vertex, and let i be the largest index such that 
v w u\ = (if no such index exists, v is the sink, and we are done). Lemmas 14.31 and 14.5 
show that after at most ( n ^~ ) — 1 iterations, we additionally have 4>( v )tt(i) = 4>( v )it(2) = 
• • • = ^(w) 7r (j_i) = + (since the next pivot step flips v w ^). At this point v is at level 1, 
since there is a unique (_): at coordinate n(i). After the flip at coordinate n(i), we re- 
peatedly apply Corollary 14.41 together with Lemma [431 to also produce l's at coordinates 
ir(i — 1), ... , n(i — 2), ... , 7r(l), at which point we have reached the sink. This takes at 
most 3(i — l)(n — l)/2 more iterations, summing up to a total of 



+ 3(i-l)(n-l)/2 < ( ) +3(n-l)(n-l)/2 = 2n 2 -(5n-3)/2. 



D 



A better bound in the case of the identity permutation can be achieved by thoroughly 
examining the run of the algorithm, as we do next. 

The identity permutation. Let us now go on to prove Theorem 14.21 For this, we first 
identify certain "milestone" vertices that are visited by Algorithm l3.1l with R = MuRTYjd, 
and then we count the number of iterations to get from one milestone to the next. Let 
us define 

v i := (0, 1, 0, 1 . . . , 0, 1, 0, 0, . . . , 0), i = 0, . . . , (n - l)/2. 

V v ' 

2i 

4.6 Lemma. // Algorithm \3.1\ with R = MuRTYid starts at the vertex v l , < i < 
(n — l)/2, it gets to the vertex v l+l in 4i + 3 iterations. 

Proof. The algorithm starts at(+:J: + :J:^+:J: ££"!), and then proceeds by the rules 
of pivoting ([6]) and (|7|) in four phases (recall that Murty^ always pivots on the left- 
most — ): 

• It pivots on the coordinates % + 1, 2i — 1, . . . , 1 and thus reaches the vertex 
( \ I \ I I | 1 + ° ::: ) after i + 1 iterations. 
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• It pivots on the coordinates 2, 4, 6, . . . , 2i, 2i + 2 and thus reaches the vertex 
(+ + + + :::+£ - + •*•'•') after * + l iterations; 

• It pivots on the coordinates 2i, 1i — 2, . . . , 2 and thus reaches the vertex 
( I |. I \ ;;; 1 + i + !!! ) after i iterations; 

• It pivots on the coordinates 1, 3, 5, ..., 2i + 1 and thus reaches the vertex 
( % \ % \ '.'.'. % \ \ \ ::. ) = v i+1 after i + 1 iterations. 

Therefore it takes Ai + 3 iterations to get from v l to v l+1 . □ 

The previous lemma allows us to count the number of iterations from = v° to 
■yl"- 1 )/ 2 The following lemma takes care of the remaining iterations. 

4.7 Lemma. If Algorithm \3J\ with R = MuRTYjd starts at the vertex v^ n ~ l >/ 2 , it gets 
to the global sink in (n + l)/2 iterations. 

Proof. The proof is similar to the proof of Lemma 14.6} except that here only the first 
phase takes place: if the algorithm starts at the vertex ( + jj_ ^. .J. ,',',' + ^ - ) , it pivots on 
coordinates n, n — 2, . . . , 1 and thus reaches the sink after (n + l)/2 iterations. □ 

Proof of Theorem \4-%\ By Lemmas 14.61 and 14.71 the overall number of iterations from 
the source to the sink is 

(n-3)/2 

v- , s n + 1 n 2 + 1 „ 

E (4* + 3) + — = ^-. D 

i=0 

5 K-matrix LCP 

In this section we examine K-LCPs, the linear complementarity problems LCP(M, q) 
where M is a K-matrix. 

5.1 Definition. A K-matrix is a P-matrix so that all off-diagonal entries are non- 
positiveo 

We introduce two simple combinatorial conditions on unique-sink orientations and 
prove that one of them implies that all directed paths from are short, and the other 
one implies that all directed paths in the orientation are short. We show that K-LCP ori- 
entations satisfy these conditions and conclude that a simple principal pivoting method 
with an arbitrary pivot rule solves a K-LCP in linearly many iterations, starting from 
any vertex. 

First, the uniform orientation is the USO in which v — > v © i if and only if V{ = (in 
other words, all edges are oriented "from to 1"). 

5.2 Definition. A unique-sink orientation <p is 2-up-uniform if whenever U = {u® I : 
I Q {hi}} is such that Uj = Uj = and u is the source of U, then the orientation of the 
subcube <f>[U] is uniform (see Figure Ej). 
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uffi j ^u © {i,j} 

A 

4>[U] i 

I 



■u®i 



Figure 3: 2-up-uniformity: The orientations of the solid edges imply the orientations of 
the dashed edges. The top row contains vertices with jth entry equal to 1 and 
the bottom row vertices with jth entry equal to 0; the left column contains 
vertices with iih. entry equal to 0, and the right column contains vertices with 
ith entry equal to 1. 

A K-USO is a unique-sink orientation that arises (as described in Section [2]) from the 
LCP(M, q) with some K-matrix M and some right-hand side q. 

5.3 Proposition. Every K-USO is 2-up-uniform. 

Proof. The essential fact we use is that a subcube of a K-USO is also a K-USO, as 
was observed already by Stickney and Watson [3U Lemma 1] . Thus it suffices to prove 
2-up-uniformity for K-LCPs of dimension 2. Hence suppose that M = (" ^) is a 2 x 2 K- 
matrix, so a, d > and b, c < 0. Then M _1 = (_f ~ ) /(ad — be) > and nonsingular. 
Now is the source of the orientation induced by the LCP(M, q) if and only if q < 0. 
Hence —M~ 1 q > 0, which proves uniformity. □ 

One particular nice property of a 2-up-uniform orientation is that all paths from the 
vertex to the sink t have the shortest possible length, equal to the Hamming distance 
|i| of t from 0. To prove this statement, we use the following lemma. 

5.4 Lemma. Let (ft be a 2-up-uniform USO and let v be a vertex such that Vi = 1 and 
v © i —> v. If Vj = and u = v © j is an out-neighbor of v, then u © i — > u. 

Proof. Suppose for the sake of contradiction that x>i = 1, Vj = and v © i — > v , and 
that u = v © j, v — > u and u — > u © i. Let f/:={«®I:K {i, j}} and consider the 
subcube 4>[U] (see Figure HI left): we observe that v © i — > u © i because otherwise the 
subcube would not contain a sink. But then it follows from 2-up-uniformity of that 
u © i — > u, a contradiction. □ 

5.5 Proposition. Let (ft be a 2-up-uniform USO and let t be its sink. Then any directed 
path from to t has length \t\. 



Another common name in the literature for a K-matrix is Minkowski matrix. 
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Figure 4: Illustrating the proof of Lemma 15.41 (left) and Lemma 15.111 (right). 



Proof. It follows from Lemma l5.4l bv induction that if u is a vertex on any directed path 
starting from 0, and if u% = 1, then u@i — » u. Hence all edges on a directed path from 
to t are oriented from to 1, and therefore the length of a directed path from to t 

is \t\. □ 

As a direct consequence of Propositions 15.31 and 15.51 we get the following theorem. 

5.6 Theorem. Every directed path from to the sink t of a K-USO has length \t\ < n. 

□ 

What is the actual strength of this theorem? We know that from any vertex of a 
unique-sink orientation there exists a path to the sink of length at most n (Lemma I3.2|) . 
It has also long been known how to find such a path from the vertex in the case of a 
K-USO [29]. The novelty is that any directed path starting from reaches the sink 
after the least possible number of iterations. Hence a simple principal pivoting method 
with an arbitrary pivot rule finds the solution to a K-LCP in at most n iterations. 

In view of this result, it is natural to ask whether it also holds if the path starts in 
some vertex different from 0. Imposing 2-up-uniformity is insufficient, as we can observe 
by looking at the Morris orientations. If we swap 0's and l's in any Morris orientation, 
we get a 2-up-uniform orientation. Hence all directed paths from to the sink are short; 
indeed, is the sink of the orientation after swapping, so all such paths have length 0. 
Nevertheless, long directed paths (and even directed cycles) do exist in this orientation. 
Therefore we introduce the following stronger combinatorial property of USOs, which 
turns out to be satisfied by K-USOs as well. 

5.7 Definition. A unique-sink orientation (j) is 2-uniform if it is 2-up-uniform and the 
orientation (j)" n >> constructed from <fi by reversing all edges is also 2-up-uniform. 

5.8 Proposition. Every K-USO is 2-uniform. 

Proof. If cj) is induced by the LCP(M, q) for some K-matrix M, then ^>(M) is induced by 
the LCP(M, —q), hence it is also a K-USO and therefore it is 2-up-uniform. □ 

Surprisingly, the simple property of being 2-uniform enforces a lot of structure on the 
whole orientation, as we show next. 
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5.9 Theorem. The length of every directed path in a 2-uniform USO is at most In. 
Before we prove the theorem, let us point out an important corollary. 

5.10 Corollary. Algorithm \3.1\ with an arbitrary pivot rule R, starting at an arbitrary 
vertex of the corresponding USO, finds the solution to an n-dimensional K-LCP in at 
most In iterations. □ 

As a consequence, we get a result for a larger class of LCPs. If M is a principal pivotal 
transform [SJ Section 2.3] of a K-matrix M', then M is a P-matrix [8, Theorem 4.1.3], 
and the USO <j> induced by the LCP(M, q) is isomorphic to the USO <fi' induced by the 
LCP(M',<j') for suitable q'. It follows that Corollary [53D] also applies to the LCP(M,q), 
even though M is not necessarily a K-matrix itself. 

Our proof of Theorem 15.91 is based on the following crucial fact, which extends 
Lemma 15.41 Informally, it asserts that once we have a (,) in some coordinate, we 
will always have a ( ). 

5.11 Lemma. Let 4> be a 2-uniform USO and let v be a vertex such that Vi = 1 and 

v © i — > v . If u = v © j is an out-neighbor of v, then u © i — > u. 

Proof. If Vj = 0, the claim follows from Lemma 15.41 So let us suppose that we have 
Vi = 1, Vj = 1 and v © i — »■ v, and that u = v © j, v — ► u and u — > u © i. Let 
U := {u © I : I C {i,j}}. We observe that i;©2—>u©i so that the subcube 4>[U] 
contains a sink. But then it follows from 2-up-uniformity of (j)(™ that v — > v © i, a 
contradiction (see Figured! right). □ 

Proof of Theorem \5.9l Let i be a fixed coordinate. Suppose in the considered directed 
path there is an edge v © i — > i> such that Uj = 1. As a consequence of Lemma 15. Ill all 
vertices uona directed path starting at v satisfy U{ = 1. It follows that in any directed 
path no more than two edges appear of the form v © i — > v for any fixed i: possibly one 
with m = and one with m = 1. This fact then implies that the length of any directed 
path is at most In. □ 

The strength of 2-uniformity can perhaps be explained by its being equivalent to what 
we call local uniformity. A USO <fi is locally up-uniform if for every U = {u © I : I C J} 

such that iij = and u —+ u @ j for all j € J, the orientation of the subcube <j>[U] is 
uniform. A USO <f> is locally uniform if both ^> and 0^"^ are locally up-uniform. 

5.12 Proposition. Let 4> be a USO. 

(i) <p is locally up-uniform if and only if it is 2-up-uniform. 
(ii) <p is locally uniform if and only if it is 2-uniform. 

Proof, (ii) is easily implied by (i); and one implication of (i) is trivial. The remaining 
implication is not difficult to prove by induction on the dimension of the considered 
subcube. □ 
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Local up-uniformity allows for a slight improvement in solving a K-LCP. When fol- 
lowing a directed path starting in 0, instead of traversing one edge at a time we can 
perform a "greedy" principal pivot: jump straight to the sink of the subcube induced 
by outgoing edges. Unlike the general P-matrix case, it cannot lead to cycling. Due to 
local up-uniformity, it will never increase the number of iterations and may sometimes 
reduce it. 

Locally uniform USOs vs. K-USOs. We have shown (Propositions 15.81 and I5.12p that 
every K-USO is 2-uniform (equivalently, locally uniform). It is natural to ask whether 
the converse is also true. 

The answer to this question is negative. It can be shown that the number of n- 
dimensional locally uniform USOs is asymptotically much larger than the number of 

/ n-l \ 

K-USOs. A lower bound of 2 V L(»^— i)/2j / on the number of locally uniform USOs can be 
derived by an adaptation of Develin's construction llOj of many orientations satisfying 
the Holt-Klee condition. An upper bound of 2°' n > on the number of P-USOs (and 
thus also K-USOs) follows from Develin's proof of his upper bound on the number of 
LP-realizable cube orientations. 

At present, however, we do not know whether locally uniform P-USOs form a proper 
superclass of the class of K-USOs. 
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